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Abstract. Let M be a complete Riemannian manifold, N 6 N and p > 1. 
We prove that almost everywhere on x = (xi,...,xjv) £ M for Lebesgue 

1 ^ 

measure in M^, the measure /J.(x) = — y * 8 Xk has a unique p-mean e p (x). 

As a consequence, if X = (Xi , . . . , Xjv) is a M^-valued random variable with 
absolutely continuous law, then almost surely fi(X(uj)) has a unique p-mean. 
In particular if (X n ) n >l i s an independent sample of an absolutely continuous 
law in M, then the process e p n (u}) = e p (Xi(uj), . . . , X n (u>)) is well-defined. 

Assume M is compact and consider a probability measure v in M. Us- 
ing partial simulated annealing, we define a continuous semimartingale which 
converges to the set of minimizers of the integral of distance at power p with 
respect to v. When the set is a singleton, it converges to the p-mean. 



1. Introduction 

Finding the mean of the median or more generaly the p-mean e p of a probability 
measure in a manifold (the point which minimizes integral with respect to this 
measure of distance at power p) has numerous applications. There is not much to 
say for the mean in M. d , almost the only case where there is a closed formula, and 
the most important case as the most useful estimator in statistics when the measure 
is uniform law on a sample. For medians in R d the situation is more complicated. 
Uniqueness holds as soon as the support of the probability measure is not carried 
by a line. The first algorithm for computing e\ is due to Weisfeld in [53]. As for the 
computation of (the center of the smallest ball containing the support of the 
measure), Badoiu and Clarkson gave a fast and simple algorithm in [5J. For many 
applications in biology, signal processing, information geometry, extension to other 
spaces is necessary. The median in Hilbert space is computed in [9]. In nonlinear 
spaces with convexity assumptions, uniqueness has been established in [18j for the 
mean, [1] for the p-mean. Many algorithms of computation now exist. As far 
as deterministic algorithms are concerned, one can cite [T5], [H], [T3], [2] for the 
mean in Riemannian manifolds, [3] for the mean in Finsler manifolds, [11] and more 
generally [25] for the median, [5] for eoo. Stochastic algorithms avoid to compute 
the gradient of the functional to minimize. They can be found in [53], [5]. For other 
functionals to minimize, see [5]. 

In this paper we investigate the case of non necessarily convex, complete Rie- 
mannian manifolds. Our first result (Theorem 12.11) concerns uniqueness of the 
p-mean of the uniform measure on a finite set {x\, . . . , x n } of points, almost every- 
where onx = (xi, . . . , x n ) for the Lebesgue measure. This generalizes Bhattacharya 
and Patangreanu result on the circle ([7], case p = 2). See also [10] for more general 
uniqueness criterions on the circle. 

l 
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For computation of the p-mean, usual deterministic algorithms are not possible 
any more, due to the fact that the functional to minimize may have many local 
minima. So restricting to symmetric spaces we use a simulated annealing method 
with a continuous stochastic process, together with an estimation of the gradient 
to minimize via a drift moving faster and faster. With this method we are able 
to define a process which converges in distribution to the p-mean for p £ [l,oo) 
(Theorem I4.3[ an d Theorem 13.21 for more general but smooth functionals). 

The main applications are in signal processing with polarimetric signal, but also 
for the group of rotations of K™, so as to determine averages on rotations. Also 
this solves many problems of optimization which may arise in economy, decision 
support, operation research. Notice that on the circle, fast computation of the mean 
has been performed in [17] ■ In fact this is a case where a closed formula can be 
found. For general case the situation is much more complicated and the convergence 
of our processes is slower and weaker. Jump processes and algorithms related to the 
continuous processes presented here will be investigated in a forthcoming paper. 

2. Uniqueness of p-means for uniform measures with finite support 

Let M be a d-dimensional complete Riemannian manifold with Riemannian dis- 
tance denoted by p. For v a probability measure on M and p > 1, we define 

H p ^ : M R+ U {+oo}, 



Either H PtV = oo or for all y £ M, H PyU (y) < oo. In the latter case we denote by 
Qp^u the set of minimizers of H pv . When Q p . v has only one element we denote 
it by e Pi „ and call it the p-mean of v. When there is no possible confusion we let 
e p = e P , v - For x = (xi, . . . ,xn) S M , we let 



Clearly H p ^ x ^ is finite. 

Theorem 2.1. Assume p > 1 or {d > 1 and N > 2}. For almost all x £ M N , 
Qp,n(x) h as a unique element e p fl{X ) 

Remark 2.2. This theorem extends Theorem 4.15 in [25] where the same result 
has been established for p = 1 and M compact. 

Proof. We begin with the case p > 1. 

Since n(x) has a finite support, we can assume that M is a compact Riemannian 
manifold. For this a smooth modification outside a large ball is sufficient. For 
instance we can choose a radius so that the boundary is smooth, double the ball 
and finally smoothen the metric locally around the place where the pasting has 
been performed. 

So in the sequel we will assume that M is compact, with diameter L. For y £ M 
we denote by S y M C T y M the set of unit tangent vectors above y. Let 



(2.3) V = {(y,n), y £ M, n = (m, . . . , n N ), nj £ S y M, j = 1, . . . N} X [0, 2L] 



Note V is a compact smooth (N + l)d-dimensional manifold with boundary. 



(2.1) 




(2.2) 




k=l 
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Define 
(2.4) 



V -> M N 



(y, n, r) n- (exp^mn), . . . , exp y (n N r N ) 



The map <fi is onto. If £ = (xi, . . . , xn) € M N , consider y £ M minimizing H p fJi ( x y 
Then among all (n, r) such that 

(2.5) 4>(y,n,r)=x 

we can choose one so that for all k = 1, . . . , N the map s t— > exp y (sn&) is a minimal 
geodesic for s £ [0, r^]- For this choice we have 



1 " 

(2-6) H pMx) (y) = -J2 r k 



k=l 



Now since y minimizes H Pjfl / x -\, from equation (|2.6p at y and variation of arc length 
formula, we have for all u £T y M 



N 



(2.7) (5>rw)<o 

and this implies 



\fc=i 



JV 



(2.8) ^^-^=0. 
So letting 



fe=i 



(2.9) W p 



fe=i J 



and 4> P = 4>\w ^ ne restriction of to Wp, </> p is onto, on M by (|2.5j) and ((27 
By Sard's theorem, the set Ci C M N of singular values of <p p has measure 0. It 

is closed since W p is compact. 
Let us prove that the set 

(2.10) C* 2 := {(xi,...,x N ) £ M N , {x 1 ,...,x N }r\Q p ^ 7^0} 

has Lebesgue-measure 0: we can assume that for i =/= j, Xi Xj since we exclude 
0-measure sets. So the elements we consider are images by (f> p of 

(2.11) W p = {(y,n,r) £ W p , n = 0, Vfc > 2 r fc > o} . 

The set W p is a submanifold of codimension 1 of W p . Now dim W p — Nd = dim M N 
so dim W'p = dimM w — 1 and its image by <f> p is of measure in M . As a 
conclusion, C2 has measure 0. 
Define 

(2.12) C 3 :={(x 1 ,...,x N )£M N , 3i £ j s.t. x % = x 3 ) 

and C = C\ U C2 U C3. The set C is closed in M N and has measure 0. Letting 

(2.13) W p = {(y,n,r)eW p , Vfc = 1, . . .N, r k e (0, 2£)} , 
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we proved that </>p|w p is onto on M N \C. Denote <j) p — 4> P \w p - Since W p has same 
dimension as M N and W p is compact, every point x of M N \C has a neighbourhood 
V x such that 0~ 1 (V X ) = U\_ x U • • • U U mx>x where the Uj jX are disjoint open subsets 
of W p and 

(2-14) 4> P \u^ ■ U jtX -> 4> P (U jtX ) 

is a diffeomorphism. Now since M N \C is second countable we can cover it by a 
countable number of such sets V x . So to prove that the p-mean is almost everywhere 
unique it is sufficient to prove it on V x . 

For x' £ V x denote x' — (x[, . . . , x' N ), and for i € {1 . . . m x }, write 

(0k„)~V) = (yi(x'),n{(x'), . . . nUx'),ri(x'), . . .,r l d (x')). 

Let i,j € {1 .. .nij} satisfy i 7^ j. If yi(x'),yj(x') G Qp^( x ') then we have 

(2-15) H P ,n(x>) Vi(x') = H P ,n{x>) 

We can assume with the same argument as for (|2.5[) and (|2.6p that the maps 
(2.16) 7 i)fc ,^ : s i-» eiq> y . ixl) (snl(x')) and j jtk , x > ■ s i-> exp y . (xl) {sn{(x')) 

are minimal geodesies respectively on [0,rjL(a;')] and [0, ri(x')]. So letting h p : 
W p -> R, (y, ra, r) H> J]fcLi r L we nave 

^h p o((t)p\ u% J- 1 (x') =H Ptll{x , ) oy i {x l ), ^h p o((t>p\ Uoa: Y 1 (x') =H Ptll{x , ) oy j (x'). 

It is sufficient to prove that for all x' € V x , 

(2-17) hp o ((kkJ- V) = h P O^kJ'V) 

implies 

(2.18) grad^, (h p o (^^J -1 ) ^ grad x , o (c^k J -1 ) . 

Indeed with (|2.18[) we will be able to deduce that the set 

(2-19) {(x' € V x , hp o O^J- 1 - h p o {tpiu^y 1 } 

has codimension > 1 in V x and this will imply that 

(2.20) {(x' £ V x , H pMxl) o y z (x') = Hp Mx , } o yj (x')} 

has codimension > 1 in V x . 

Let us prove ([2718]) , For fc = 1 , . . . , iV let 

KOO = -Ti.fc.x'K^')) and TO i(- T ') = -ii.M'W <>'))■ 
These unit vectors satisfy 

ex Px' k ( r U x ') m k( x ')) = Vi( x ') and exP^W^')^^')) = 2/i OO • 

Then noting that (/i p o (^pkJ -1 ) (a:') = Z)fc=l( r *) P ( a; *) we S et 

da:' (0pk,J _1 ) (') 

/ N \ 



= (-p^(riT-H x >i( x '),T x , yi (.)\ 

\ fc=1 It t t\M 

- 1 y i las') 1VA 

-pdirUx'^mKx'), . . . , (rj,(a/)) p - 1 «*jvOO) - -> T 
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Due to the fact that (yi(x'), n l {x'), r l (x')) e W p , the first term in the right vanishes. 
So 

(2.21) 

gracL. {hp o (^I^J- 1 ) = -p ((rUx'^mlix'), {r^f^m^x')) 

and similarly 
(2.22) 

grad,, (hpo^lu.j- 1 ) = ^((ri^OF-'mi^),---,^^'))^ 1 ^^)) ■ 

Since yi(x') ^ Vj(x') we have (r\(x'),m\(x')) ^ (r^(x'), m\{x')), so (rj(x')) p_1 mj(a;') 7^ 
(r^(a;')) p_1 TO^(x'), from which we conclude that 

grad x , (hp o {(j) p \u^y l ) ^ grad a , (h p o (^1^. J" 1 ) . 

This achieves the proof for the case p > 1 . 

Let us now consider the case p = 1. The result is due to Yang in [2B], we give 
the proof here for completeness. 

The main difference is that the subset of M N of points x = (x±, . . . , xn) so that 
x i G Qi,fj,(x) f° r some i has positive measure. 

First consider the open subset U of M N of points x such that for alH = 1, . . . , N, 
Xi Qi,n(x)- 

Consider the closed subset Co of M N of points (x\, . . . ,xn) — <fi(y,n } r), with 
(y, n,r) g V such that for all j, k = 1, . . . N, rij — in^. Since d > 1 and N > 2 
this subset has Lebesgue measure 0. 

Replacing by J7 and C by Co U C, the argument is similar until (|2.17p . But 
now we will be able to prove that (|2.17l) implies (|2.18l) only in some neighbourhoods 
V XyX > to be precised later, of x' £ V x such that the geodesies 

s i-> exp y .( xl) (sni(x')) and s h-> exp Vj ^ n (sn> k (x')) 

are minimal respectively on [0, rlfa/)] and [0, ri(x')]. But this will be sufficient 
since every compact subset of V x can be covered by a finite number of these neigh- 
bourhoods ~V XiX t . 

Making the above assumption on x', the proof is similar until (I2.2ip and (|2.22l) . 
Then we have 

(2.23) grad x , (hi ° folk,.)" 1 ) = - (m\(x'), . . .,m? N (x')) 
and 

(2.24) grad,, (hi ° (faK.)' 1 ) = - (m{(x'), . . . , m? N (x')) . 
Assume 

grad^, (hi o (0ik,J _1 ) = grad^, (hi o J" 1 ) . 

Then for all k = 1,...,N, m k (x') = m k (x'). In particular for k = 1 this im- 
plies (possibly by exchanging i and j) that yi(x') lies in the minimizing geodesic 
from x[ to yj(x'). Now since x' ^ Co there exists fc € {1,. . .N} such that xj. ^ 
{exp v i i x '){ sn i( x ')), s G [ — 2L,2L]}. On the other hand since m l k (x') = m k (x'), 
yj(x') (or yi(x')) lies on the minimizing geodesic from rrj, to yi(x') (or yj(x')). 
As a consequence there are two minimizing geodesies from yi(x') to yj(x'). But 
this is impossible since the geodesic from x\ to yj(x') is minimizing, contains 
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Ui(x') and x[ ^ Ui{x') by the fact that we have supposed that x[ ^ Qi^( X ') an d 
yi(x') € Qi,n( X ')- So 

grad^, (hi o (0i|£/ M ) _1 ) + grad^, (/n o (^ilt/j.J -1 ) , 

and by continuity this is true in a neighbourhood V XlX > of a;'. 

Now we consider the case where x[ £ Qi )jU (x') an d ^ Qi,^(V)- We follow the 
same lines as in the previous part with the difference that now yi(x') — x[ and for 
the definition of Ui }X W\ is replaced by 

Wi = {(y,n,r)eV, n = 0}. 

The definition of Uj. x remains unchanged. By [25] Theorem 1 

1 N 



fc=2 



< M ;v(z')(K}) 



which gives 
(2.25) 



A' 



fc=2 



< 1. 



Since d > 1 and N > 2, the submanifolds of 14 images of 

AT 



(y,n,r) e t/i^, 



E 

fc=2 



= f 



and 



AT ~| 

(y,n,r) € C/^, X! ™ fc = f 

fc=2 J 



by ^1 have measure 0, so we can exclude them. On the subset 

N 

(y,n,r) e U itX , < 



< 1 



J2 nk 

k=2 

the function h\ is smooth and on its image by 4>\, 

(2.26) grad,, (hi o (^k.J" 1 ) - - (0, m\[x'), m%(x')) . 
Again 

(2.27) grad,, (hi ° (^|^ J" 1 ) = - (mj(x'), . . . , m? N (x')) . 

They are not equal, and this achieves the proof for this case by the same argument 
as before. 

Finally we consider the case where x[,x 2 £ Qip(x') with x[ — yi{x') and x' 2 = 
yj(x'). We follow the same line as in the previous case, but now for the definition 
of Uj tX , W\ is replaced by 

W{ = {(y,n,r)£V, r 2 = 0}. 
Again we can exclude the submanifolds of V x images of 



(y,n,r) € U ja 



ke{l,...,N},k=£2 



= 1 
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and 



by and work on 



(y,n,r) €U jjX , n k = 

fce{i,...,A r },fc#2 



0i S (V,n,r) e C/ ix , < 




On this set hi o (<fti\ui x ) 1 and hi o (f^il^ !C ) 1 are smooth and 

(2.28) grad,, (hi o (fcl^ J" 1 ) - - (0, m\[x'), m%(x')) . 

(2.29) gvad x ,(hio^ 1 \ U] J- 1 )=-(m[( x , ) 1 0,mi(x'),...,m%(x')) 
They are not equal, and this achieves the proof. 



□ 



Corollary 2.3. Let p £ [Loo) and X — (Xi,...,Xn) a random variable with 
values in M , which has an absolutely continuous law. Then almost- surely /i(X(w)) 
has a unique p-mean e p (X(uj)). 

Corollary 2.4. Let p € [l,oo) and (X n ) n >i a sequence of i.i.d. M-valued random 
variables with absolutely continuous laws. Then the process of empirical p-means 



(e p , n (uj) := e p (X 1 (uj), . . . ,X n (u))) ) 

\ / n>l 



is well-defined. 



Remark 2.5. For p = 2 and M a circle, it has been proved in [7] that the assump- 
tion can be weakened: the same result holds if the law has no atom. 

We believe that it would be interesting to study the behaviour of the process 
( e f),n)n>i m many situations. For instance when the law of Xi is uniform on a 
compact symmetric space (even the case of the circle is highly non trivial) one 
would observe a recurrent but irregular and slower and slower process. Again 
on a compact symmetric space, when the law v of Xi has a finite number of p- 
means due to a finite group of symmetries, one would observe an almost stationary 
behaviour, and at increasingly spaced times jumps between smaller and smaller 
neighbourhoods of the p-means of v. 



3. Finding the minimizers of some integrated functionals with 

simulated annealing 

Let M be a compact Riemannian manifold. For simplicity and without loss of 
generality we assume that M has Lebesgue volume 1 . On M consider a probability 
law v with a density with respect to Lebesgue measure, also denoted by v. Assume 
we are given a continuous function k : M x M — > R + , where n(8,y) is interpreted 
as some kind of cost for going from 6 to y. Assume furthermore that for all y E M 
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the function 9 i— > n(9,y) is smooth and that its first and second derivative in 9 are 
uniformly bounded in (9,y). Consider on M the functional 

U : M -> K+ 

(3 ' 1} 0H- / <9,y) V {d,y) 

Denote by .M the set of minimizers of U. The aim of this section is to find a 
continuous semimartingale which converges in law to M.. Also we try to avoid 
using the gradient of U, which in many cases is difficult or impossible to compute. 

For this we will use a sequence (Pk)k>o of independent random variables with 
law v, a Poisson process N t on N with intensity 7 t _1 where 

(3.2) lt = (l + t)-\ 
Define 

(3.3) c{U) =2 sup ( inf e(<j>)) , 

denoting the set of continuous paths [0, f] — > M and for <p € %,y; the elevation 
e{4>) being defined as 

(3.4) e(4>) = sup U{<j>{t)) - U{9) - U(y) + inf U(z). 

o<t<i Z ^ M 

Let 

(3.5) ft = i]n(l + i), 

the constant k satisfying k > c(U). 

We assume that (N t )t>o is independent of the sequence (Pk)k>o- We let (B t )t>o 
be a Brownian motion with values in R r for some r € N, independent of (N t )t>o 
and (Pfc)/t>o, and a a smooth section of TM ® (R r )*: for all € M, cr(0) is a linear 
map IT -> T e M. We assume that for all e M, we have a(9)a(9)* = id Tf>M - We 
fix #o G M and let t be the solution started at 9 of the Ito equation 

(3.6) d@ t =a(Q t )dB t - p t gr&d @t K(;Y t )dt with Y t = P Nt . 

Recall that if -P(Ot) : T So M — >• T Qt M is the parallel transport map along (9 t ), 
then 

(3.7) de t = p(e t )d ^ p(e s )- x o de s 

Also define 9° the solution started at 9o of the Ito equation 

(3.8) d9° = a(9°) dB t - ft Qf^ grad e o «(-, y) v{y)dy^j dt. 
Note p.8p rewrites as 

(3.9) d<d° t = er(9") dB t - ft grad e? U dt, 

so that the same equation with fixed /3 instead of ft has an invariant law with 
density 



(3.10) uptf) = ±re-^ e \ with Z = [ 



M 
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The process 0® is an inhomogeneous diffusion with generator 



(3.11) 



L° t (0) = -A(8)-f3 t gr a d g U. 



Denote by m t {9) the density of t . 

The process (8<, Yf) is Markovian with generator L t given by 

(3.12) 

Ltf(9, v) = - Pt grad e K (;y)\ /(•, y) + 7^ J (M z) - f(0, v)) v(dz) 

= L ltt f(;y)(0) + L 2 , t f(O,-)(y). 

We know that for all neighbourhood Af of M, J^j- fip(0) d0 converges to 1 as 
(3 00. So to prove that Jj^m t (0)d0 converges to 1 it is sufficient to prove the 
following proposition: 

Proposition 3.1. The entropy 

m t {9) 



(3.13) 

converges to as t — > 00 
Proof. Let us compute 
(3.14) 

dm t {9) 



Jt := I In 

M 



mt(6) d9 



dJt_ 
dt 



A I 



dt 



d9 - 



M 



dln(i f3t (9) 
dt 



m t (0) d9 



In 



(m t {0)\ dm t {0) 



AI KVPttf) 



dt 



■ d0. 



Since for all t mt(0) is a probability density, the first term in the right vanishes. So 
we get 



(3.15) 



dJt 
dt 



= 2!3' t / U(0)(m t (9)-^M)d9 + 



A I 



A I 



In 



( m t (0) 



m t (9)d0 



where the last term comes from Dynkin formula. For the first term in the right we 
have using (|3.5[) 

(3-16) 2(3' t J U(0)(m t (0) - ^M) d9 < 4||[/|U#| < ^M^. 

Now we split the second term in the right of (|3.15[) into 
mt{0) 



(3.17) 



A I 



hi 



= / L» 

M 



In 



We have 



L° t 



A I 



In 



( m t {9) 
WW 

(m t {0) 



m t {9) d9 



m t {9)d9+ / Rt(0,y) 

M 



In 



m t {9) 



m t {0) d0. 



(3.18) 



L 



A I 



WW 

m t {0) 



A I 



m t {9) 



m t (0)d0 
tip t (0) d0~ 

MM® 



A I 



Vln 



f m t (8) 



< -2c 2 ((3 t V l)- p cxp (-c{U)fo) J t 
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for some C2 > and integer p > by logarithmic Sobolev inequality ( [15j and 
[16], for more details see [2D])- Note we used again Dynkin formula to prove the 
vanishing of the first term in the right of the second line. 
As for the second term we have 



R t (e,y) 



A I 



hi 



Pt{ din 
din 



m t (fl) 
m t (0) 



m t (0) d9 
, grad e 



f m t {6) 

m [ J^(9)(dJ^(9) 1 R t (9))m t (9)d9 

Jm V m « \ V ma / 



rad e «;(•, z) v{dz) ) mt(9, y) d9 dy 



grade «(•>£/) ("^(j/! ) ~ ) w t (0) ^ dy 



with 



Rt(0) 



Rt(9,y) 



A I 



hi 



<2& 



-/ 




Jm 




inequality 


m t (6) \ 










2 











grade K ('> y){ m t{y\0) - v{y)) dy. 



m t {9) d9 



1/2 



M 



<pt / Ilii t (0)ll' ! m t (0)d0 + 



M/9t 



|i? t (0)|| 2 m t (0)<20 



1/2 



Replacing 2c2 by C2 in (|3.18[) we can after summing get rid of the second term in 
the right. Defining 



(3.19) 

let us now bound 

\R t {9)\\ 2 m t (9) d0 = 



K= sup || gradg «;(•, t/) || 

0,yeM 



M 



A I 



gr&d g K(-,y){m t {y\9) ~ v{y))dy 



Al 



m t {9) d9 



< 



K 



I M 

< 32K 2 
= 32K 2 I t 



\m t (y\9) - v{y)\ dy 



A I 



m t (0)d0 



Al 



where we have defined 
(3.20) I t 



In 



dy. 



m x m \ v {y) J 

We also used classical bound of total variation by entropy ([IB]): 



\m t (y\9)-v(y)\dy<4V2 



hi 



A I 



rrif 



(y\Q) \ 



m t (y\9) dy 



1/2 
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At this stage we proved that 



(3.21 



dJ t < 4||k| 



c 2 (/3 t V l)- p ex P (-c{U)/3 t ) J t + /3 2 32K 2 I t . 



dt ~ k(l + t) 

The next step is to find a suitable bound for As before 
dl t 



(3.22) 



dt 



Li 



M x M 



in (^m 

. V nv) J 



'MxM 

We begin with the first term 



/ (L 2 ,t + L ht ) 

J MxM 



m t (y, 9) dOdy 

m t (y,9) dOdy. 



V nv) J 



MxM 

It' 

-1 
It 



. V nv) J 



MxM J M 



In 



m t {y,6) dOdy 
( m t (z\6)\ T _ ( m t (y\6) \ 



In 



Ky) / 



v(dz)mt{y, 9) dOdy 



In 



mxm V v {v) J 



By Jensen inequality we have 



In 



MxM 



mi 



(y\Q) \ 



/m VJm ny) 



( m t (y 
m t (y\9) 



v{y)m t {6) dydO 



m \Jm V nv) J J 



v{y)dy)m t {9) d9 



M 



ln(l)m t (0) d9 = 



Consequently 



La 



MxM 



<-7f 



In 



Ky) / 



m t (y,0) dOdy 



In 



which rewrites as 
(3.23) 



£2 



MxM 



mxm \ v \y) J 
m t (y\6)\ 



d6 



In 



Ky) / 



m t (y,e)dOdy<-jf 1 I t . 



Let us now consider the second term in the right of (|3.22|) . Since 



In 



it rewrites as 



(3.24) 



MxM 



hi 



m t 



(Q\y) \ 



mt{6) ) 



m t (y,6) dOdy 



L 1>t [ln(m t (%)) - Hm t (6))] m t (y, 0) d9dy. 



MxM 
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But 
(3.25) 



MxM 
1 



L M ln(m t (Q\y))m t (y,Q) d8dy 



MxM 



MxM 



Aln(m t (%)) - p t / (detnmt(-\y),gr&dgK(-,y))m t {y,d)dddy 



M 



Vy/mt(9\y) deu{dy) 



- Pt (d e m t (-\y),gr&d e n(-,y))v(y)d8dy. 

J MxM 

Let us bound the absolute value of the last term: 



-fit I (d g m t (-\y) 7 gY&d e K(-,y))v(y)d9dy 

I MxM 



2/3* / {deV m t(-\y), gradg V m t(0\y)v(dy) dOdy 

MxM 



(3.26) 



MxM 



VvW(%) y/mt{e\y) V {y)dSdy 



< 



This yields 
(3.27) 



MxM 



MxM 



-fiK 2 m t (0\y)+2 VvW%) )v{y)d6dy 



MxM 



L lit ]n{mt(e\y))m t (y,6)d0dy < -ffK 2 



We also have to bound the last term in p.24[) : 

/m f (fl) 



(3.28) 



-L M ln(m t (0)) = -L ljt ln 



WW 



l m Mma(0)). 



We already know that 
(3.29) 



£ llt In f ^) m t (0) d9<-^+ 4|| K |U^. 

For the second term we have 

-L M (0)) =2L ltt U(6) 
(3.30) - ftAC/(0) + /3 t 2 (d£7, grad e «(•, I/)) 

< K'ifit V 1)A 

with 
(3.31) 



K'= sup |A eK (-,2/)|+^ 2 . 

8 ,y£M 



Finally we obtain 
(3.32) 
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together with (pHTjt : 

(3 ' 33) l£ - WTT) ~ ° 2{Pt V l)" P «p(-^)A) J * + 2/? t 2 32X 2 / t . 

At this stage we can use the end of the proof of theorem 1 in 21 to obtain that 
under assumptions (|3.5[) and (|3.3[) then 

(3.34) lim J t = 

(notice that in Section |4] we will prove this in a more general context). □ 
Theorem 3.2. Assume 

(3.35) ft = -hn(l + t), and 7t = (l+t)- 1 , 

where k > c(U), (c(U) defined in (|3.3|) ). Then for any neighbourhood Af of M, 

(3.36) lim P [Q t e Af] = I. 

t—>oc 

Proof. We use Proposition 13.11 together with the fact that 

\\mt-»p t \\ < 4v/2Jt 

and 

lim Hfrl/f) = 1. 

t— >oo 

a 

4. Application to location of p-means in symmetric spaces 

In this section we assume that M is a compact symmetric space endowed with the 
canonical Riemannian metric of volume 1. Denote by p the Riemannian distance 
in M, D its diameter. We fix p > 1 and consider a probability measure v on M. 
We aim to find at least one element of Q PtlJ by using the result of the previous 
section. In particular if v has a unique p-mean e p , then we will be able to construct 
a process which converges in probability to e p as t — > oo. 

Denote by p(s, x, y) the heat kernel on M, and for s > let v s be the probability 
measure with density 

(4.1) v s (y) = p(s,y,z)i/(dz), 

J M 

and let 
(4.2) 
and 
(4.3) 



■ s : M x M -> M + 

(6,y) i y / p{s,0,z)p p {z,y)dz 1 

J M 

[/.„.,„ : M — >• R-l 



8 ^ K Sl (9,y)u S2 (y)dy. 
Jm 

Also let U = H pM . Clearly u Sl and n S2 satisfy the assumption of the previous 
section. Moreover, denoting by A4 SliS2 the set of minimizers of U Sl . S2 then as 
s\,S2 — > we have A4 Sl , S2 — > Qp,u is the sense that for any neighbourhood Af of 
Qp.v, we have A4 Sl ,s 2 C A/" for all si,sa sufficiently small. This is due to the fact 
that as si,S2 — > 0, U Sl ,s 2 (8) — > U{9) uniformly in 9. 
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Lemma 4.1. For all si,S2 > we have 

(4.4) U auS2 (6) = Uo,s 1+ sM= [ P (0,y)v Sl+S2 (y)dy. 

J hi 

Proof. Fix 6,y € M, let m be the middle point of a minimal geodesic from 9 to y 
and i m the symmetry centered at m. We have 



p(s 1 ,9,z)p p (z,y)dz = / p(s 1 ,i m (6'),i m (z))p p (i m (z),i m (y))dz 

J M 

= / p{s\,i m {9),z')p p {z' ,i m (y))dz' 

J M 

= [ P (s 1 ,y,z')p p (z' ) 6)dz' 

J M 

pP(0,z')p(s 1 ,z / ,y)dz' 



' M 

where we first used the invariance by isometry of the heat kernel and then did the 
change of variable z' = i m (z) in the integral and finally used the symmetry of the 
heat kernel. To finish the proof we are left to use the convolution property of the 
heat semigroup. □ 

Corollary 4.2. Defining 

(4.5) K"= sup || grad e p(.,y)||, 

6,y£M 

we have for all s±, s 2 > 0, 9,y £ M , 

(4.6) || grade k s1 (-,j/)|| <pD p ~ 1 K"=:K and || grad e U sliSa \\ < K. 

With all these properties we would like to find s\(t) \ and s 2 (t) \ such 
that the process ©t started at 9$ and solution to 

(4.7) d&t = a(e t )dB t -p t gr a d etKsi(t) (-X 2it) )dt 

converges in law to e p , (N t ,Y t S2 ^) being a Poisson point process in [0,oo) x M 
with intensity r y{t)~ 1 v S2 (t){y) dt dy, independent of (B t ). This is the object of the 
next theorem in which we will take 

s 1 {t) = s 2 {t)=s t = {\n{l+t)y 1 . 

So define 0° the solution started at 9o of the Ito equation 

(4.8) d<3° t = a(e° t ) dB t - ft (J grad e o K St (; y) u„ t (y)dy^j dt. 

Notice that using Lemma |4~TI (|4.8p rewrites as 

(4.9) dQ ( l = (7(9°) dB t - ft grad e? U 2st dt, 

where U 2st '■— Uo^ 2st , so that the same equation with fixed (ft s) instead of (ft, St) 
has an invariant law with density 

(4.10) ^ s {6) = J-e-W-W W ith Z Pt .= [ e- 2 ^'Ue>. 
The process 0^ is an inhomogeneous diffusion with generator 

(4.11) L° t (6) = ±A(0) -Pt grad 9 U 2st . 
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Denote by mt(0) the density of Gt. 

Let Yt := Y t St . The process (&t,Yi) is Markovian with generator Lt given by 

(4.12) 

Ltf(0, y) = Qa(0) - ft grad e n St (•, y)) /(•, y) + 7^ (f(9, z) - f(9, y)) v St (dz) 

= L lit f(;y)(9) + L 2it f(6,-)(y), 

We know that for all neighbourhood M of Q Pll /, jj^ np tS (9) d6* converges to 1 as 
j3 — > oo, uniformly in s sufhciently small. Again define 

J f := / In ^ 



(4.13) J t := / In m t (0) d0. 

Theorem 4.3. Assume 

(4.14) ft = 1 ln(l + i), 7t = (1 + t)- 1 , «i(t) = a 2 (t) = *(*) = (ln(l + t))" 1 - 

where k > c(U), (c(U) defined in (|3.3|) ). T/ien /or an?/ neighbourhood Af of Q p ^, 
the process Qt defined in equation (14. 7f> satisfies 

(4.15) 

Proof. We use Proposition 14.41 below together with the fact that 

IK-MA, St || < 4^/2Jt 

and 



Km P[6 t € A/] = 1. 

t— >-oo 



lim fi0 t)St (N) = 1. 

t— s-oo 



□ 



Proposition 4.4. T/ie entropy 



(4.16) 



m t (9)d9 



converges to as i — > oo . 
Proof. Let us compute as before 

(4.17) ^ = -[ d t H^ uSt (e))mt(0))d9+ f L t 
al Jm Jm 

For the first term in the right we have using (|4.10[) 
(4.18) 



In 



mt(g) 



m t (9) dO. 



-2[3' t U 2st - 2ft 



MxM 



2s' t d s lnp(2s t , 6, z)p(2s u 9, z)p p (z, y) v(dy)dz 



+ 2(3' t / U 2st (9>0 t ,s t (9')d9' 
Jm 

+ 2ft f ( f 2s' t d s \np(2s t , 6', z)p(2s u 9', z)pP{z, y) dzv{dy) ) np tlSi (9') d6' . 
It is known that there exists Cq > such that Vs £ (0, 1] 



(4.19) 



\d s \np(s,9,z)\<^-, 
s l 
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see e.g. [14] and [22] where bounds of the type |Vg In p(s. 9. z) \ < — and |V \np(s, 9, 
C 2 

— are given. Here we use 



\d a \np(a,e,z)\ 



1 



A g p(s,e,z) 



p(s,9,z) 



dimAf r, , , , , , 
< —5— (\V 2 g ]np(s,d,z)\ + |V e lnp( S ,#,z)| 2 ) 



So (j418l) and (|4J9j) yield 



(4.20) ' </;p( 4 ^ + ^r^ 



which implies 
(4.21) 
with 

(4.22) C = DP(4 + C ). 

Evaluating with (|4. 14[) and integrating on M we get 



(4.23) 



5 t ln(^ t , St (0))m t (0) d9 



M 



< 



c 



(l + t)k 



(l + ln(l + <)). 



Now we split the second term in the right of (|4.17j) into 
(4.24) 



Li 



M 



In 



M 



In 



m t {6) 
( m t (9) 



m t {9,y) dOdy 



m t {6)d0+ I R t (9,y) 

M 



In 



f m t (9) 



We have as for (|3TT51) 



(4.25) 



M 



In 



m t (0) 



m t (9)d9 = -2 / 

J M 



m t (6) 



m t (0, y) e?0dy. 



< -2c 2 (/3 t V l)- p exp(-c(C/ 2st )/3t) J t 



for some c 2 > and integer p > by logarithmic Sobolev inequality ( |20) ) . 

The computation for the second term is similar to the one after (|3.18|) and we 
get 



/ Rt(B,y) 



In 



( mt(9) 



m t (9,y)d9dy 



2/3* 



M V m t 



nil 



(9){dJ -{9),R t (6))m t {9)d6 



with 



Rt{9) 



;rad e K St (-,y){m t (y\6) - v„ t (y))dy, 



M 



and again 



Rt(9,y) 



M 



In 



m t {9) 



m t (6) dO 



\ \\R t (9)\\ 2 m t {9)d9- 



M 



M 



VJ^-[<f) 



MusM do. 



Replacing 2c 2 by c 2 in (I4.25P we can after summing get rid of the second term 
the right. 
Here again 



/ \\R t (9)\\ 2 m t (9)d9 <32K 2 I t 

J M 



where we have defined 
(4.26) It = 



( m t {y\6) \ 
In — m t {y,6)dy. 



At this stage we proved that 
(4.27) 



djf C , , . . 
< — (l + ln(l + f)) 



dt ~ {l + t)k 

- c 2 {f3 t V l)- p exp (-c{U 2st )Pt) Jt + /3 2 32K 2 I t . 

The next step is to find a suitable bound for As before 



(4.28) 



and 



(4.29) 



Now 



-4- = - I d t \n{v st (y))m t (9, y) d9dy 
dt J M 



Lt 



MxM 



In 



(y\9) \ 



Vs t {y) J 



m t (y,9) dOdy 



d t ki(v St (y))m t (9, y) dQdy 



M 



< 



s't / d 8 \np(s t ,y, z)p(s t ,y, z) v(dz) 

J M 

Wt\C _ c 



(4.30) 



MxM 



s 2 



hi 



l + t 



(y\9) \ 



v St {y) J 
/ (La lt +Li 

J MxM 

We begin with the first term: 



m t (y,0) d9dy 



In 



m, 



(y\9) \ 



s t (y) J 



m t (y,9) d9dy. 



L, 



MxM 



It 



In 



mt 



(y\9) \ 



v St {y) ) 



m t {y,Q) dOdy 



In 



fmtm\ {UsM „ mtm)mt[e)de 



mxm V v s t {y) J 
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and estimate it as for (|3.23[) : 



(4.31) 



MxM 



In 



v s t {y) ) 



m t (y,8)d6dy < -7^% 



For the second term in the right of (|4.28j) we need to introduce the density ft (y) of 



hi 



( m t (y\0) \ ( m t {9\y) \ 
\vsMJ \m t {6) ) 



hi 



My) 
v s t (y) 



it rewrites as 



(4.32) 



Li, 



MxM 



hi 



nil 



(8\y) \ 



mt(6) J 



m t (y,0) d6dy 



L 1>t [ln(m t (%)) - ln(m { (0))] m t (y, 9) dOdy. 



MxM 



Similarly to (15^5]) 
(4.33) 

Li t t \n(m t {Q\y))m t (y, 9) dOdy 



MxM 
1 

2 



MxM 



Aln(m t (9\y)) - ft / (dg lnm t (-|y), grad e K St (■, y)) m t (y, 9) dOdy 



M 



= -2 



MxM 



VvW%) Mft(v)dy 



(d e m t (-\y), grad e K s t (-,y)) f t (y)d6dy. 



'MxM 

For the absolute value of the last term: 



(4.34) 



-ft 

1 



(d e m t (-\y), grad e «(-, y)) v(y) f t (y) dOdy 



MxM 



< ^/3 2 K 2 + 2 



MxM 



VVm*(%) My)d8dy. 



We get as in (|3T27|) 

(4.35) 



'MxM 

Then we bound the last term in (I4.32[) : 



L 1>t \n{m t {e\y))mt{y,e)d6dy < -j3 2 K 2 



(4.36) 



-Lx tt \n{m t {0)) = -L M ln 



( m t {d) 



We already know by (jiTf]) and (pT2"5|) that 

m t (6) 



Li :t ln(iJ,/3 uH (e)). 



(4.37) - / L M ln 
Jm 



mt(0)d9<-^ 



C 



dt (1 + t)k 



(l + ln(l+f)). 



For the second term we have 

L u \n{ti Pt , St (6)) = -2ptLi, t U 2a M 
(4.38) = -ft AC/ 2st (6) + 2ft 2 (dU 2st , grad e K St (• , j/)> 

< X'(ft V l)/3 tSt " 2 
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for some K' > 0, where we used 

AU 2s = f (A e lnp(2 S , 6, y) + \\V e lnp(2s, 6, y)\\ 2 ) p(2s, 6, y)ff(y, z) u(dz) 
Jai 

and standard bound for the first and second derivatives of the heat kernel ([14] and 

[22])- 

Finally we obtain 

(4.39) + fJcKl + *) V k)(Hl + t)f - f - (1 + ,)/« 

together with P~2"T)l : 

(4.40) 
j r 

< (1 + ln(l + t)) - c 2 {p t V !)-*> exp (-c(f/ 2st )/3 t ) J t + 2f3 2 32K 2 I t 



which rewrites as 



(4-41) -± < h (ln(l + t)f - - (1 + t)I t 

and 

(4-42) ^ < Cl + (ln(l + t)fl^ c 2 (ln(l + <))^(1 + t)"^ J, 

for some constants ci, fci > 0, as soon as f > 2. At this stage we can use a similar 
computation to the end of the proof of theorem 1 in |21j to obtain that under 
assumptions (|4. 14[) and (|3.3[) then 

(4.43) lim J t = 0. 

t— J-oo 

However we will do the calculation for completeness, and because there are some 
small differences. Recall U s —> U uniformly as s — >• 0. Moreover 2st — > as t — >• oo, 
so we get 

limsupc(C/ 2;5t ) < c(U). 
As a consequence, for t sufficiently large we have 

(4.44) <l- £ 

K 

for some e > 0. Let 

ci(ln(l + t)) 2 



(4.45) £ t = 



l+t + ci(ln(l + t)f - c 2 (ln(l + t))-P(l + *)-(*-«) 



where e > is defined in (|4.44j) . It is easily checked that for t sufficiently large it 
is positive and decreasing, and that it converges to as t — > oo. Define 

(4.46) K t = J t + i t If 

We will prove that K t — > as t — > oo and from this we will get (|4.43|1 . 
for t sufficiently large, 
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and this yields with (|4.41|) and (|4.42|) 

- £ t ci(ln(l + t)fl t - (1 - ^)c 2 (ln(l + t))- p (l + t)- Cj ^J t 
+ £ t k 1 (\n{l + t)) 4 -{l + t)£ t I t . 

Replacing ci(ln(l + 1)) 2 at the end of the first line by 

£t(l + t + Cl (ln(l + t)f - c 2 (ln(l + t))- p (l + 

by the help of (|4.45l) we obtain 

dK t In(l + t) „ , ,, „, , <v2s t ) 
< cx \ - t > - c 2 (, t (\n{l + t))-P(l + t) — I t 

- (1 - 4)c 2 (ln(l + t))-P(l + t^-^Jt + * t fci(ln(l + i)) 4 
and this yields using -(1 + t) « < -(1 + i)~ (1_e) : 

(4.48) ^ < A t - fl t tf, 
with 

(4.49) A t = ci ^ 1 ^ + £th ln(l + 1)) 4 
and 

(4.50) B t = (1 - ^)c 2 (ln(l + + 1)-* 1 -*). 
A sufficient condition for K t to converge to as t — > oo is 

/OO 
B t eft = +00 

and 

(4.52) lim — = 0. 

Condition (|4.51l) clearly is realized. As for condition (|4.52l) we easily see that 





ln(l+t) 
C1 1+i 







(l~£t)c 2 (\n(l + t))-P(l + t)-(^) 

and also 

^fc!(ln(l + f)) 4 

(l-£ t )c 2 (ln(l + i))-J>(l + *)-(!-*) 
from the fact that 

< c(ln(l+t)) 2 
* " 1 + i 

for some c > 0. □ 
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